livio montanaro 18/08/2022

formula completa semiarchi placidus senza latitudine pianeta
complete formula semiarchs placidus without planet latitude

non traduco in inglese, le formule sono intuibili.
ob 	= obliquità 23.438130

dati forniti dal software
	lob/lo[0][a]	= longitudine pianeta promissore
	lat		= latitudine 
	dp      	= declinazione

dh      = distanza oraria
ad      = differenza ascensionale
a       = significatore (0 = sole 1 = luna da tabella effemeridi)
aa      = promissore

aspectangle[aaa] = aaa = tutti gli aspetti es 0,60,90,120 ecc

distance=command pd corrisponde all'angolo minore fra due punti
normalize= command pd riduce entro 360°

:psig0;`calcolo significatore xlat0 parte 1
u=1;
lob= lo[0][a];
los =lob;
ED = sin(los)*sin(ob);
ED = asin(ED);`distanza equinoziale Gansten
da[u] = tan(lat)*tan(ED);da[u] = asin(da[u]);
Sad[u] = normalize(90 + da[u]);semiarco diurno
San[u] = normalize(90 - da[u]);semiarco notturno
dp = decp[a];
lp = 0;@cra;	calcolo ar senza latitudine
dm_mc[u] = normalize(ar-armc);if (dm_mc[u] > 180) dm_mc[u] -= 360;distanza meridiana mc
dm_ic[u] = normalize(ar-aric);if (dm_ic[u] > 180) dm_ic[u] -= 360);distanza meridiana ic
return;

:ppro0;`promissore xlat0 parte 2
ora si calcola l'angolo del promissore alla congiunzione o agli aspetti col significatore
partendo dalla longitudine del promissore,aggiungendo o sottraendo l'angolo

lob = normalize(lo[0][aa]+/-aspectangle[aaa]);` nuova longitudine del promissore

'la direzione è sempre calcolata con latitudine 0 (lp = 0) proiettata sull' equatore

u=0;
lop = lob;
ED = sin(lop)*sin(ob);
ED = asin(ED);`gansten
da[u] = tan(lat)*tan(ED);da[u] = asin(da[u]);
Sad[u] = normalize(90 + da[u]);
San[u] = normalize(90 - da[u]);
dp = decp[aa];
lp = 0;@cra;
dm_mc[u] = normalize(ar-armc);if (dm_mc[u] > 180) dm_mc[u] -= 360;
dm_ic[u] = normalize(ar-aric);if (dm_ic[u] > 180) dm_ic[u] -= 360);
if (normalize(-lo[0][14] + los) > 180)
{`se sopra orizzonte
dh = dm_mc[1] / Sad[1]; 
delta = normalize(dm_mc[0] - (dh * Sad[0]));
}
else
{`se sotto orizzonte
dh = dm_ic[1] / San[1];
delta = normalize(dm_ic[0] - (dh * San[0]));
}
delta = angolo direzione finale
si consiglia di usare 1 grado = 1 anno
return;